In-gap impurity states as the hallmark of the Quantum Spin Hall phase 



J. W. Gonzalez 1 , J. Fernandez-Rossier 1 ' 2 . 
(1) International Iberian Nanotechnology Laboratory - INL, 
Av. Mestre Jose Veiga, ^115-SSO Braga, Portugal 
(2) Departamento de Fisica Aplicada, Universidad de Alicante, San Vicente del Raspeig, 

(Dated: March 4, 2013) 



Spain 



We study the different response to an impurity of the two topologically different phases shown 
by a two dimensional insulator with time reversal symmetry, namely, the Quantum Spin Hall and 
the normal phase. We consider the case of graphene as a toy model that features the two phases 
driven, respectively, by intrinsic spin-orbit coupling and inversion symmetry breaking. We find that 
strictly normalizable in-gap impurity states only occur in the Quantum Spin Hall phase and carry 
dissipationless current whose quirality is determined by the spin and pseudospin of the residing 
electron. Our results imply that topological order can be unveiled by local probes of defect states. 

PACS numbers: 



The intrinsic properties of an electron gas are revealed 
in the way it reacts to the presence of a localized pertur- 
bation. In metals, the period of Friedel oscillations pro- 
vides information about the Fermi surface [1]. In semi- 
conductors, the binding energy of shallow acceptors and 
donor states depends on the effective mass and dielectric 
constant of the host material [2]. In superconductors, 
both the impurity induced modulations of the density 
of states Q and the presence of zero energy or mid-gap 
states [4], reveal the symmetry of their order parameter. 
Here we address the fundamental question of whether 
two dimensional topological insulators react to a local- 
ized perturbation in a way different from conventional 
insulators. 

Topological insulators [5] have a bulk bandgap like a 
normal insulator but have protected conducting states on 
their edges and surfaces. This bulk-boundary correspon- 
dence has been rigorously established [6|, l7j in boundaries 
that preserve the translational invariance in at least one 
dimension. In two dimensional topological insulators, 
edge state are expected to present quantized conductance 
[8MlO| which can be used to unveil the existence of bulk 
topological order [11]. More recently, the bulk-boundary 
correspondence has been extended to the case of topolog- 
ical defects that lead to Hamiltonians H(k, r) that vary 
slowly with adiabatic parameters r surrounding the de- 



fect [12]. 



Here we explore the electronic structure of a two di- 
mensional insulator, that can be either in the Quantum 
Spin Hall (QSH) or in the normal phase, in the neighbor- 
hood of an isoelectronic impurity that creates a repulsive 
short range potential. We find that only in the QSH 
phase in-gap states appear and have exotic electronic 
properties: they carry non-dissipative spin current. Our 
finding provides an alternative way to detect topological 
order, using local probes sensitive to the density of states 
in the neighborhood of the impurities. 

We use gapped graphene as a toy model for two dimen- 
sional topological insulators [8|, |9[ . The electronic prop- 



erties of graphene are intimately related to the structure 
of the honeycomb lattice, formed by two interpenetrat- 
ing triangular sublattices, A and £?, related by inversion 
symmetry, which define a pseudospin degree of freedom 
that we denote with the operator r z . There are two ways 
to open a gap in graphene preserving the size of its two 
atom minimal unit cell. Both lead to interesting elec- 
tronic phases. A conventional gap opens in graphene 
when a sublattice symmetry breaking potential, ^r Zl 
is included in the Hamiltonian. This gap entails pecu- 
liar electronic properties: the two valleys carry orbital 
currents of opposite sign [13]. When the gap is opened 
by intrinsic spin orbit coupling, as described with the 
second-neighbour spin dependent hopping proposed by 
Kane-Mele, graphene is in the QSH phase. 

The Kane-Mele @,0| tight-binding Hamiltonian Hq de- 
scribes electrons in a hexagonal lattice, with first neigh- 
bor hopping t, spin-dependent second neighbor hopping 
t so and the ^-t z term. This model commutes with S z , 
the spin projection perpendicular to the graphene plane, 
and each spin sector is identical to the Haldane model Q 
for spinless fermions. We consider the effect of a substitu- 
tional isoelectronic impurity, described with a single-site 
potential, in an otherwise boundless and perfect gapped 
system: 



H = H + VoJ2\ 0s )(° s l 



(1) 



The strength of the spin-independent impurity potential 
is Vo, and acts only on the atom of the A sublattice. 

For each spin channel s, the crystal Hamiltonian Ho 
can be written in the form of a 2x2 matrix in the sublat- 
tice basis: 

W k = t [/(k)r+ + /* (k)r-] + ^t z + st so g(k)T z , (2) 

where r = (r x ,T y ,T z ) are the Pauli matrices in the sub- 
lattice space, r ± = r x ±ir y , s = ±1 are the eigenvalues 
of the spin operator S z . /(k) and g(k) are the usual 



functions that sum the Bloch phase over the nearest and 
path-dependent next-nearest neighbors [9j. 

The Bloch Hamiltonian is identical to that of a 
(pseudo)spin r in an effective field: T~L^ = • r, where 

Ao+t so g(k)~ 



t/(k)+tr(k),t/(k)-tr(k), 



The energy bands for the model are given by e I/jS (k) = 
±e^, where v = ± labels the two bands (per spin chan- 



nel) and 



|^| . In the top panel of Fig. [T]we plot 



them for a given value of s, along the line that joins the 
Dirac points K and K' , all of them with the same values 
of t and t so = O.lt and different values of Aq. Time rever- 
sal invariance ensures that, for the opposite spin orienta- 
tions, we have e I/)S (k) = k). The K and if' points 
define the so called valley index a = ±1. These points are 
special because the in-plane components of the effective 
field vanish, making the energy splitting between the 
bands minimal, 
with: 



At K and K' we have hi 



(0,0, 7«0 



4(Ao 



scrAi), 



(3) 



where scr can take only two values scr — ±1 and Ai = 
6-s/3£ so - Thus, at the Dirac point the wave functions have 
a well defined sublattice (r z ) character. The effective field 
at the Dirac points, ±7 SC r, defines the top of the valence 
band and bottom of the conduction band at each valley. 
We plot them in the low panel of Fig. [TJ 

In the transition from the QSH phase, with 7_ < 
to the normal phase, with 7_ > 0, the system closes the 
gap (7_ = 0), reflecting the impossibility of deforming 
adiabatically one phase into the other [8, 9]. In the QSH 
phase, the energy splitting is finite at both valleys, the 
orientation of the effective field h^. is opposite at K and 
K' for a given s. In this phase, the model presents topo- 
logically protected edge states inside the gap, making the 
edge metallic. At the normal phase, with 7_ > the ori- 
entation of the effective field is now the same at both 
valleys and a gap opens in the edge states. 

We now explore the possible appearance of in-gap 
states in the two insulating phases of the model in the 
presence of the impurity potential, as described by Eq. 
(PQ) . This is different from previous works that have stud- 
ied the influence of impurities on the conducting surface 

We address 
17|. We de- 



MUM- 



states of 3D topological insulators 
our problem using the T-matrix formalism 
fine the Green function operators, Q{z) = 
and Go(z) = (zI — Ho) -1 where z is a complex number 
and I is the unit matrix in the Hilbert space of the infi- 
nite lattice. A closed expression for the Green function 
Go(z) can be readily obtained in terms of the eigenstates 
of H^.. For the on-site spin-independent potential, the 
complete Green function is related to Go(z) by: 

V 



G(z) = G°(z) + G°(z)\0)- 



VoGUz) 



(o\e°(z), (4) 
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FIG. 1: (Color online) Top panel: energy bands, for the s — + 
spin channel, for t so = O.lt, and 4 values of Ao. Bottom panel: 
evolution of bands at K and K' and in-gap impurity energy 
Eb(Ao) calculated both with the lattice model (black) and 
the continuum model (red). 



where Go o(z) is the unperturbed Green function diagonal 
matrix element in the atomic representation and is the 
atom at which the perturbation is located. 

The appearance of bound states with energy would 
be given by the existence of poles of G in the band-gap 
defined by H^.. Thus, we have to solve the equation 

1 

Vo' 

In order to have a closed expression for Go o(^>) we use 
the Lehman representation and project over the site rep- 
resentation. Any given site in the lattice i can be iden- 
tified by its unit cell R{ and the sublattice T{ = A,B. 
The unperturbed Green function is written in terms of 
the eigenstates of the crystal Hamiltonian 



(5) 



\v, k,s) 



R,T 



where N is the number of unit cells of the crystal, 
and IAI k T are the components of the eigenstates of 
the Bloch Hamiltonian in Eq. (|2]). If we express 
the effective field in spherical coordinates as, h^. = 
e£. (sin#£ cos sin 0£ sm <£k> cos ^k)' ^ ne corresponding 
wave functions for the v = ± bands read: 



sin I - I e 



^ U S -, K b = -™s[-j. 



K,KB = sin - , (6) 
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where we have omitted the subscripts from 6 and cp for 
the sake of clarity. The unperturbed Green function ma- 
trix can be written as: 

1 _ pk.(Ri-Rj) 

= nE z _ £k „ feJ ^, (7) 

Using this expression in combination with Eq. J6j), we 
can recast the Eq. JHJ) as 



1 

Vo 



2N 



E 



1 



nt - z 



1 



E h + E b - el 



(8) 



where ng. = ^k/ e k- Importantly, the Green function 
(?oo ( z ) is expressed as a sum over the whole Brillouin 
zone of a function that depends on the projection of the 
effective field vector and, as such, it contains information 
of the topology of the Bloch states, which is a necessary 
condition to expect a relation between the solutions of 
Eq. J5J) and the topological order in the system. 

The Eq. JHJ) can be analytically solved in the strong 
coupling limit V^ 1 = 0, in two cases. For Ao = 0, the 
repulsive potential yields a mid-gap state E b = whose 
properties we discuss below. For Ai = 0, there is a 
solution with Eb = — Ao/2, in agreement with a general 
result [18] . However, this solution is degenerate with the 
top of the valence band and it is a resonance rather than 
an actual in-gap state. The interpolation between these 
two limits is obtained by the numerical solution of Eq. 
(JHJ) and is shown in Fig. [I] for Vo = I0 6 t. Bound states 
are only found when Vo > t and, interestingly, only when 
Ao < Ai, i.e., in the topological insulator phase. This is 
the main result of the paper: we find that a local impurity 
can bind an in-gap state only in the QSH phase. 

In order to obtain some analytical insight of the one to 
one relation between the existence of in-gap states and 
the topological phase, we have worked out Eq. J7J) in the 
continuum limit, in which only states close to the two 
Dirac points are included. Their crystal Hamiltonian is 
then given by 



Ho(k : a,s) = Hv F (k x r x + iak y r y ) + —r z 



Ax 



sar z . (9) 



For this model it is possible to obtain a closed analytical 
expression for Eq. JHJ): 



1 

Vo' 



2ir(hv F y 



^2(<y S T-E b ) log 



1 



(10) 

where e c is the cutoff energy [19]. The sum over sr re- 
flects, for a given spin orientation, the contributions com- 
ing from both valleys. In the strong coupling limit, V Q ~ X 
is negligible and the existence of an in-gap solution of 
Eq. JTOj) requires that the right hand side sum vanishes. 
Since 7 + — E b is always positive, E b must satisfy 7_ < E bl 
which, by definition of in-gap state, is only satisfied in the 




FIG. 2: (Color online) Local density of states as a bar in 
a exponential scale. The in-set shows the impurity-induced 
current between links for the state at Eb. The impurity site 
is shown in yellow. 



QSH phase and not in the normal phase (see lower panel 
of Fig. [TJ). Thus, the continuum model also has the one 
on one relation between topological order and the emer- 
gence of in-gap impurity states. The numerical solution 
of Eq. [TOj for V ~ l = is shown in Fig. [TJ 

Now we address the electronic properties of the topo- 
logical in-gap impurity states. The Green's function for- 
malism provides a closed expression for the in-gap wave 
function in terms of their binding energy E^, obtained 
from Eq. J7J), and the unperturbed Green function [l7| : 



\4> b ) = 



dQUE b ) 
dE 



Y,G%{E b )\i). (11) 



We consider the case A = and V^ 1 = 10 -6 t. Expect- 
edly, the in-gap wave function is localized around the im- 
purity site, as shown in Fig. [2j Interestingly, the in-gap 
states is not fully sublattice polarized, in contrast with 
zero energy states in bipartite Hamiltonians [io|. Thus, 
we have ((j>b\(7 z\<j>b) — 0.56, significantly below 1. We also 
find, analytically, that the wave function in the impu- 
rity sublattice is purely imaginary, while for the other 
sublattice is purely real. In the reciprocal space it is 
also true that the wave function has unequal weight on 
both valleys, but is not fully valley polarized. Both the 
incomplete sublattice and valley polarizations are quan- 
titatively different from spin-filter edge states in zigzag 
ribbons for the same model, which are fully sublattice 
and valley polarized. 

In contrast, the most salient feature of the edge states 
is also present for the topological in-gap impurity states: 
they carry current, which is quite unusual for a local- 
ized state. The current operator is defined at the bonds 
of the tight-binding Hamiltonian imposing the continuity 
equation. For a given pair of sites n and m in the lat- 
tice, the current operator reads [2l[ J nm = j nm |n)(m| — 
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jmn\m)(nl where j mn et nm /ih and t nrn is the Hamilto- 
nian matrix element (n\Ho\m). Thus, the current carried 
by bound states reads 

I^ d = ^Im[t nm <P* b (m)Mn)}. (12) 

The map of the current is shown, for a given spin, in 
the inset of Fig. [2j for t so = O.lt, A = 0, V = 10~ 6 t. 
Although the in-gap state is more localized in the B sub- 
lattice, opposite to the impurity site, the current is larger 
in the AB bonds than in the BB bonds because t = 10t so . 
For the opposite spin, the current flow changes sign so 
that the net current is zero but the spin-current is not. 
We have also verified that, when the impurity site is in 
the other sublattice, the current flow is inverted. Thus, 
the in-gap states have dissipationless spin-currents whose 
chirality is determined by the sublattice at which the im- 
purity resides. 




20 40 60 80 100 

FIG. 3: (Color online) Bound, band and total current in an 
AB bond , as a function of Vb. Upper inset: density of states 
on B sublattice. Lower inset: evolution of Eb(Vo) for Ao = 0. 

The observable current probed experimentally would 
be given by the contribution of all occupied states, which 
includes both the in-gap state and the valence band, 
whose density of states are shown in the upper inset of 
Fig. 03 The contribution to the current from the band 
states reads 

7- 

^ =^ J Im [Gnm(E)j mn - g mn (E)j nm ] dE. (13) 

— oo 

where Q is the full Green function whose closed form is 
given in Eq. (j4]). Thus, there are two contributions to the 
current, one given by the band states and the other given 
by the in-gap states. We plot both them in Fig. [3j for 
A = 0, as a function of the impurity strength, Vq. They 
have opposite signs and, in the strong Vq limit, cancel 



each other. The evolution of £^(Vb) is shown in the lower 
inset. In that limit the impurity bound state becomes a 
mid-gap state Eb = and the spectrum recovers electron- 
hole symmetry, for which ground state currents are not 
possible [22| . The cancellation of the edge current due to 
the contribution of the bulk states takes also place in the 
case of ribbons. However, a finite non-zero spin-current 
is obtained in our case for a wide range of Vb close to 
what it is expected for absorption of atomic hydrogen in 
graphene. 

It must be stressed that, in contrast to in-gap states 
in the domain wall of polyacetylene, the impurity in-gap 
states do not give rise to charge fractionalization. This 
can be understood as follows. Together with the in-gap 
solution of Eq. (|5]) where E^ ~ V^ 1 , there is always a 
second solution with energy E^ ~ Vq. Thus, there are 
two states outside of the bands, and not only one as in 
the case of Su-Schrieffer-Heeger solitons [23], so that both 
bands lose a complete state. 

In conclusion, we have studied the electronic properties 
of the reaction of a two dimensional insulator, as given 
by the the Kane-Mele model, to a strong on-site impurity 
potential. We have shown that in-gap impurity states 
appear in the case topologically non-trivial or Quantum 
Spin Hall phase jH, El- Vacancy induced in-gap states 
have been also predicted for a two-dimensional topolog- 
ical insulator described with a modified Dirac equation 
[24} . Here we have shown the one-on-one relation be- 
tween the topological order and the existence of the in- 
gap states in the context of the Kane-Mele lattice model. 
We have found that these topological in-gap impurity 
states carry a net spin-current which, due to the single 
occupancy of the state, is upgraded to a net current or 
orbital magnetization. Therefore, we propose that the 
topological nature of this type of insulators can be stab- 
lished by local probes of the electronic properties of these 
defects, as opposed to the highly non-local probes re- 
quired to test quantized transport of the edge states. Fur- 
ther work should determine if our results can be extended 
to higher dimensions and to other models for topological 
insulators in two dimensions. 
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